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SN 2002BU - ANOTHER SN 2008S-LIKE TRANSIENT 
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Draft version August 15, 2012 

ABSTRACT 

We observed SN 2002bu in the near-IR with the Hubble Space Telescope, the mid-IR with the Spitzer Space 
Telescope and in X-rays with Swift 10 years after the explosion. If the faint L# ~ 10 2 Lq HST near-IR source at 
the transient position is the near-IR counterpart of SN 2002bu, then the source has dramatically faded between 
2004 and 2012, from L ~ 10 6 °L to L ~ 10 4 ' 5 Lq. It is still heavily obscured, Ty ~ 5 in graphitic dust models, 
with almost all the energy radiated in the mid-IR. The radius of the dust emission is increasing as R oc t ' A 
and the optical depth is dropping as Ty oc f" L3±0 4 . The evolution expected for an expanding shell of material, 
Ty oc f" 2 , is ruled out at approximately 2er while the Ty oc f" 8 to f" 1 optical depth scaling for a shock passing 
through a pre-existing wind is consistent with the data. If the near-IR source is a chance superposition, the 
present day source can be moderately more luminous, significantly more obscured and evolving more slowly. 
While we failed to detect X-ray emission, the X-ray flux limits are consistent with the present day emissions 
being powered by an expanding shock wave. SN 2002bu is clearly a member of the SN 2008S class of tran- 
sients, but continued monitoring of the evolution of the spectral energy distribution is needed to conclusively 
determine the nature of the transient. 

Subject headings: stars: evolution - stars: supergiants - supernovae:individual (SN 2002bu) 



1. INTRODUCTION 

Sup ernova (SN) 2002bu was discovered on 2002 March 
28 by iPuckett & Gauthierl (120021) in the gal axy NGC 4242 , 
with a fairly low peak magnitude My ~ -15 (Hornoch 2002). 
The last pre-discovery observations on 2001 February 21 and 
March 14 placed a unfiltered magnitude limit at the location 
of the SN of 20.5 mag. A low resolution spectrum taken by 
lAyani et af] (120021) on April 1 showed a flat continuum and 
strong, narrow Balmer emission lines (FWHM ~ 1 100 km/s) 
that led to a Type Iln classification of the SN. The early (~80 
day s post peak) light c urve presented by iFolev et all d2007l) 
and ISmith et all d201 ll) showed a plateau and a slower rise 
and decline as compared to other Type II SNe, with the color 
b ecoming redder with tim e . 

[Thompson et al. (2009) proposed that SN 2002bu may be 
a member of a new clas s of stellar transients , with the proto- 
types being SN 2008S dArbour & Bolesll2008l) and th e 2008 
optical transient (OT) in NGC 300 dMonardl [20081) . The 
class is characterized by a transient-progenitor pair, where 
the transient is 2 - 3 mag fainter than a regular core collapse 
SN (ccSN), with narrow emission lines and evidence of in- 
ternal extinction in the spectrum, while the progenitor is a 
dust enshrouded, low luminosity (^5x\Q a Lq) s tar with lit- 
tle variability in th e years before the outburst (Prieto et al. 
120081 lPrietol l2008b . In mid-IR color-magnitude diagrams 
(CMD), the progenitors occupy the extreme end of the asymp- 
totic giant branch (AGB) sequ ence and the obscuring dust is 
graph i tic rather than sili cate (Prieto et al. 2009, Wesso n et alj 
2010). Thomps on et al.l d2009t) further showed that these tran- 
sients are relatively common among ccSNe (^20%), but their 
progenitors are extremely rare (<^ 10~ 4 of massive evolved 
stars), which implies that many massive stars go through this 
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dust obscured phase shortly (< 10 4 yrs) before the explosion. 
The rarity of the progenitor star s was furthe r confi rmed in the 
survey of additional galaxies bv lKhan et alj d2010l) . 

Since the main characteristic of the new class is a dust en- 
shrouded progenitor, it is impossible to unambiguously clas- 
sify SN 2002bu as a member, as there are n o pre-explosion 
IR ob servations of the region. However, Thompso n et alj 
(2009) found a bright, red mid-IR source at the location of 
SN 2002bu in Spitzer Space Telescope (SST) data taken 2 
years after the explosion, which indicates dust formation. 
This, together with the transient characteristics (low luminos- 
ity, dust and narrow emission lines visible in the spectrum), 
mak e SN 2002bu a like ly member of the new class. In addi- 
tion, [Smith et al] d201 lb analyzed 5 spectra taken between 1 1 
and 8 1 days after discovery, and they show that the spectrum 
of SN 2002bu becomes redder with time and evolves from a 
spectrum resembling an LBV in outburst at early times, to one 
more similar to the spectra of SN 2008S and the NGC 300 OT 
8 1 days later. The evolution of the Ha line from a Lorentzian 
profile to an asymmetric, blueshifte d Gaussian profile, is sug- 
gestive of newly formed dust. In iKochanek et alj d2012l) we 
argued that the early light curve also indicates dust formation, 
albeit with some amb i guities . 

In IKochanek et alj d2012l) we analyzed archival Hubble 
Space Telescope (HST) and SST data from approximately 
two years after the tr ansient peak. As previously reported by 
Thom pson et al.l (12009ft. SN 2002bu was a luminous mid-IR 
source in 2004, but IKochanek etal J d2012l) also found that it 
was invisible in archival HST data from 2005 in the BVRI 
bands to limits of ~ 25 mag. The spectral energy distribution 
(SED) in 2004 was well fit by surrounding ar,= 20000 K, 
L* ~ 1O 5 9 L0 source with an optically thick (Ty — 30) shell of 
dust located at a radius of approximately R ~ 10 15 8 cm. 

Broadly speaking, there are three possible explanations of 
the observations, as we discussed in our mid- IR survey of the 
"supernova impostors") Kochaneket al. 2012). The first, "tra- 
ditional", view of these events is that a shell of material is 
ejected during the optical transient and forms dust once the 
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FIG. 1. — The top two panels show the [3.6] /im SST observations of 
SN 2002bu from April 2004 (left) and January 2012 (right), the middle two 
panels show the corresponding [4.5] fim SST data, and the bottom two panels 
show the Fl 10W (J) and F160W (H) HST observations from February 2012 
of the region marked with a square on the SST images. The four top panels 
are 12"0 X 12"0. The two bottom panels are 3"0 X 3"0 (about 80 X 80 
parsecs) and the radius of the circle is ~0"3, which is 3 times the uncertainty 
in the astrometry. 



ejected material becomes cool enough. As the shell expands, 
its optical depth drops as r oc 1 ft 2 and the characteristic tem- 
perature drops as T oc L(f) 1 / 4 ?" 1 / 2 because the shell expands 
with the ejecta velocity R ~ v s t. The dust radius at the time 
of the previous SST observations of SN 2002bu was consis- 
tent with a shell expand i ng at the velocity of v s = 893 km/s 
adopted by ISmith et ail (1201 ll) . The second possibility is 
that the optical transient is a signal that the star is entering 
a high mass loss phase with a dense wind that forms dust and 
obscures the source. While the wind is steady, the optical 
depth is roughly constant with a dust temperature close to the 
dust destruction temperature Tj ~ 1500 K unless the wind be- 
comes optically thick in the mid-IR. When the high mass loss 
pha se ends, the evolution q uickly resembles the first scenario 
(see IKochanek et al.ll2012l). T he third scenario is the one in- 
troduced by IKochanek! d201 11) to explain SN 2008 S and the 
NGC 300 OT. Here the progenitors are already shrouded by a 
very dense wind when an explosive transient occurs that de- 
stroys most of the dust to leave the transient little obscured at 
peak. The wind is so dense, however, that the dust reforms 
and re-obscures the transient. The present day luminosity is 
a combination of a surviving star (if any, nothing in the data 
requires one) and the luminosity generated by the shock prop- 
agating through the wind. At later times, the optical depth 
outside the shock is dropping as r oc 1 jt and, once the optical 
depth is low enough, the X-rays produced in the shock should 




FIG. 2.— A 2"0 X 2?'0 region of the F160W HST images with a cross 
marking the estimated position of SN 20()2bu. The arms of the cross span 2<r 
of the astrometric uncertainties. The circles mark the three sources detected 
by DOLPHOT within the 3<r error in the estimated SN location with A being 
the closest and C the most distant. 

be observable. 

These scenarios make different predictions for the time 
evolution of the transient, so monitoring the evolution of 
SN 2002bu as a function of wavelength should reveal the na- 
ture of the event. Here we report near-IR HST, mid-IR SST 
and Swift X-ray observations of SN 2002bu taken roughly 10 
years after the transient peak and 8 years after the the last HST 
and SST observations. Section 2 presents the new observa- 
tional data. We discuss the results in Section 3, and consider 
their broader implications in Section 4. 

2. OBSERVATIONS AND DATA ANALYSIS 

We observed SN 2002bu with SST in June 2011 and Jan- 
uary 2012, in both the [3.6] and [4.5] ^m bands (program ID 
80015) with exposure times of 240 sec (8 dithered 30 sec ex- 
posures) for both bands. We also observed it with HST in 
February 2012 using the WFC3/IR camera and the F110W 
(J) and F160W (H) filters (proposal ID 12450) with exposure 
times of 2x700 s for each band. Figure Q] shows the changes 
in the [3.6] and [4.5] images between 2004 and 2012 as well 
as a zoomed in view of the region in the / and //-band HST 
images. The source has become significantly fainter in the 
mid-IR, more so at [3.6] than at [4.5], but any near-IR coun- 
terpart must still be very faint. Figure[2]shows a close-up view 
of the region in the F160W image. We identified 12 reason- 
ably isolated stars in the HST images that could be matched 
to the SST images using the IRAF immatch package. Using 
the 2004 3.6/im image as the reference frame, we obtain the 
estimated SN position in the HST data shown by the circle 
in Figure Q] and a cross in Figure [2] The uncertainty in the 
position is approximately 1.5 HST pixels or 0."09. 

The Spitzer fluxes were measured with aperture photometry 
( apphot package in IRAF ), following the procedure described 
in IKochanek et alj ([2012). We measure the flux of the source 
using a range of source and sky apertures (with appropriate 
aperture corrections) and then combine the individual results 
taking into account both the formal statistical and the system- 
atic uncertainties implied by the scatter in the fluxes found for 
different aperture combinations. We analyzed the HST im- 
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ages using the DOLPHOT dDolphinll20"00l) photometry pack- 
age with its standard WFC3 parameter files. DOLPHOT iden- 
tifies three sources within a radius three times larger than our 
estimated astrometric uncertainties. We label these sources 
A, B and C in order of their distance from the estimated po- 
sition (see Figure The photometric results are presented 
in Table[T]toget her with the previous H ST and SST flux mea- 
surements from Kochane k et al.l (120121) . 

We adop t a dis tance to NGC 4242 of 5.8 Mpc based on 
iTully et aTl (12009|) and a foreground Galactic extinction of 
E(B-V) = 0.01 (Schlegel et alJJ1998l) . We note, however, that 
the distance to NGC 4242 is uncertain. Th ere is a second 
Tully- Fisher distan ce estimate of 10.4 M pc by Springob et al. 
(2007) (erratum in Springob et al. 2009), and if we associate 
NGC 4242 with the group containi ng NGC 4258 it lies a t 
the intermediate distance of 7.2 Mpc (Herrnstein et al. 1999). 
Our qualitative conclusions are unaffected by these uncertain- 
ties. Quantitatively, luminosities increase proportional to d 2 , 
dust radii increase as d, and the velocities implied by the dust 
radii increase as d. The actual photometric models are es- 
sentially distance independent other than these scalings. We 
comment on the effects of distance changes as necessary. 

Figure [3] shows near-IR H/J-H CMDs constructed from 
the DOLPHOT catalogs for a large (2f3x2.'0) and a small 
3. "5 (100 pc) radius region around the SN. We required 
the signal-to-noise ratio in both filters to be greater than 4 
for a detection (see iDalcanton et aTl 1201 ll) . and a sharpness 
parameter sharpness 2 < 0.1 to exclude non-stellar sources. 
The CMDs are corrected for Gal actic extinction {E(B - V) = 
0.012 m ag. ISchlegel et all 119981) and assume a distance of 
5.8 Mpc (Tully et al. 2009). We mark the locations of sources 
A, B and C in both CMD s. The right panel in Figure [3] also 
shows the Padova (Mari go et al J 120081) isochrones for 10 75 , 
10 8 and 10 83 years, that have ZAMS masses corresponding to 
their end points of 9, 5 and 4M Q , respectively. The isochrones 
are moderately bluer than the stars, which could be evidence 
for E(B-V) ~ 0.4 mag (for 5.8 Mpc, 0.2 mag for 10.4 Mpc) 
of local extinction. Unfortunately, the optical detection limits 
are not strong enough to constrain these possibilities. For all 
these distances and extinctions, the local stellar population is 
only consistent with the older isochrones, ages closer to 10 8 
years and a maximum mass of 5M Q than 10 7 5 years and a 
maximum mass of 9M Q . Since adding these extinctions have 
little effect on our near/mid-IR models, we add no additional 
extinction beyond Galactic to the models in §3. 

We also observed SN2002bu with Swift dGehrels et all 
2004) between April 26-30 201 1, for a net exposure of 15 ks. 
We reprocessed the XRT data using the xrtpipeline 
tool provided by the Swift team and reprojected the ob- 
servations into a single image. We chose the source re- 
gion to be a circle centered on the SN with a radius of 10 
pixels (23".6) and a nearby background region without any 
sources. We did not detect the source either in the full 
(0.2-10 keV) or the 0.2-0.5 keV, 0.5-2 keV, and 2-10 keV 
bands. We obtained la limits of < 5.4, < 1.3, < 3.8, and 
< 3.4 x 10~ 4 cnt s" 1 , respectively, for these bands, where 
we corrected the s mall aperture used in the analysis based 
on the Swift PSF dMoretti et al .1120051) . Assuming thermal 
bremsstrahlung emission with Tx = 0.8 keV and Galactic ab- 
sorption of NHCaiactic = 1 • 17 x 10 20 cm" 2 (Di ckey & Lockmanl 
1990), we obtained a flux limit of < 0.9 x 10" 14 erg cm" 2 s" 1 
in the 0.5-2 keV band that roughly corresponds to a luminos- 
ity limit of Lx < 10 4 L©. More generally, we model the emis- 



sion as thermal brehmsstrahlung emission at a temperature Tx 
obscured by an additional column density Nh associated with 
the source. We use the PIMMS model for the Swift/XRT to 
estimate the absorbed soft, medium and hard-band count rates 
and constrain the unabsorbed total luminosity to be consistent 
with the observed upper limits for the three energy bands. 

3. RESULTS AND DISCUSSION 

We used DUSTY dlvezic & Elitzuri 119971; llvezic et al l 
119991) to fit the SEDs of the sources detected in HST near- 
IR data. Figure |4] shows the probability averaged SEDs of 
SN 2002bu in May 2004 and January 2012 (for sources A, B 
and C, treating near-IR HST measurements of source A ei- 
ther as detections or upper li mits). The best fit m odel to the 
2004/2005 data is taken from lKochanek et all (120 12b and cor- 
responds to a = 20000 K, L» ~ 10 s 9 Lq star surrounded 
by an optically thick (tv ~ 30) shell of dust expanding at 
^900 km/s. Both graphitic and silicate models give similarly 
good fits to the data. 

We also embedded DUSTY in a Markov Chain Monte 
Carlo engine to model the SEDs, varying the dust tempera- 
ture, Td, optical depth, Ty, and stellar temperature, T*, with 
either a fixed 2: 1 ratio between the inner and outer radii of the 
shell or allowing the ratio to vary between 1.1 and 10. We 
included a weak prior on the stellar temperature, log^T* = 
4.0 ± 0.3, (and restricted its range to 3000 K < < 30000 K) 
and on the expansion velocity implied by the inner radius of 
the shell, log 10 (v s /km/s) = log 10 (893) ± 0.30 = 2.95 ± 0.30. 
For the 2012 epoch we ran models considering star A as either 
a detection or an upper limit. The results allowing for varia- 
tions in the shell thickness were little different from those with 
a fixed thickness and no particular thickness was preferred, so 
we only report the results for the fixed thickness in Table[2] 

The results for the 2004 epoch a re the same as in our ear- 
lier models from Kochanek et al. (2012). The source must 
be quite luminous, ~ 1O 5,92±OO2 L0, but with an indeter- 
minate source temperature because of the heavy obscuration. 
The optical depth in graphitic models is log 10 Ty — 1.6 ±0.2 
(scattering plus absorption), and it is moderately higher in 
the silicate models because of the higher s cattering opaci- 
ties of silicate dusts (see the discussion in Kochanek et al. 
d2012l) ). Significantly higher optical depths begin to have sig- 
nificant opacity even for the shorter wavelength IRAC bands, 
inconsistent with the shape of the SED. The dust tempera- 
tures at the inner edge of Tj ~ 1 100 ± 170 K are roughly 
in the range expected for newly forming dust. The dust ra- 
dius is estimated to be log 1() (R, n /cm) ~ 15.77 ± 0.12, imply- 
ing a velocity of log 1() (v s /km/s) ~ 2.95 ± 0. 12 that is consis- 
ten t with the v e locity prior taken from the line widths cited 
by iSmith et alj d201 lb . Note, however, the uncertainties are 
much smaller than those of the prior, so the data are in fact 
determining this radius and velocity rather than the prior. De- 
spite the reasonable coverage of the mid-IR SED, graphitic 
and silicate dusts fit the data equally well. As discussed ear- 
lier, the primary effect of distance uncertainties is simply to 
rescale the luminosities, velocities and distances. In particu- 
lar, adopting the distance to NGC 4258 (7.2 Mpc) raises the 
vel ocity from v it ~ 900 km /s to v s ~ 1100 km/s, while using 
the ISpringob et al. (2009) distance of 10.4 Mpc raises it to 
v s ~ 1500 km/s. One reason we adopted the smaller distance 
is that the larger distances begin to require significantly larger 
mean expansion velocities than im plied by spe ctroscopic ob- 
servations during the transient (see Smith et al. 2 0TH) . 

The interpretation of the SED in 2012 depends critically 
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FIG. 3. — Near-IR CMDs from the HST F110W and F160W observations, corrected for Galactic extinction of 0.012 and assuming a distance of 5.8 Mpc 
iTully et al. 2009). The left panel shows a large region (2?3x2'0) while the right panel shows a ~3"5 (100 pc) radius around SN 2002bu. Th e locations of 
the three objects in the SN 2002bu region (A, B, and C, see Fig. [2) are marked with labeled squares. The curves on the right panel show Padova I Ma rigo et alj 
2008) isochrones for 10 7,5 , 10 s and 10 8 ' 3 years, alternating solid and dashed lines, which have end points corresponding to ZAMS masses of 9, 5 and 4Mq 
respectively. The vertical arrows shows how much the points would shift if instead we adopted the distance of 10.4 Mpc. The diagonal arrows show how the 
points would shift after correcting for E(B — V) = 0.4 mag of additional extinction while keeping the distance of 5.8 Mpc. 



on whether any source in the HST images corresponds to 
SN 2002bu. If source A (or B/C) is a detection, then the 
present day source must be a relatively hot, > 15000 K, 
star. For graphitic dust, the star is relatively low luminosity 
L, = 10 4 5±0 2 L Q , with moderate log 10 r y ~ 0.70 ±0.25 ob- 
scuration, at a relatively large radius log 10 (f?,„/cm) ~ 16.18 ± 
0.25 that corresponds to a lower average expansion velocity 
log 10 (v s /km/s) ~ 2.72 ±0.25 but is uncertain enough to be 
consistent with no change in velocity. The dust temperature 
at the inner edge is much cooler, ~ 500 K. In the sili- 
cate models the source is more luminous, = 1O 4,95±O18 L0 
with greater optical depth, log 10 Ty — 1.21 ±0.12, located at 
a smaller radius log 10 (7?,„/cm) ~ 15.98 ± 0.17, that now re- 
quires a slower expansion rate log 1() (v s /km/s) ~ 2.52 ± 0.17. 
The dust temperature at the inner edge is moderately warmer, 
T& ~ 600 K. There is no basis for choosing between the 
two dust types. If the association with source A is simply 
a coincidence and we instead treat the HST fluxes as upper 
limits, then solutions with cool stars are allowed, but they 
require increasingly high optical depths for lower tempera- 
tures, ry ~ 100 for the coolest (T* = 2500 K) models. This 
is essentially the limit of the fully obscured progen i tors of 
SN 2008S and the NGC 300-OT dPrieto et all I2008L Iprietol 
200811 iThompson et"alll2009l lKochanekll20TTlK 

The current position of the source in the mid-IR CMD (see 
Figure |5J, close to the tip of the AGB sequence, is similar 
to that of the SN 2008S or NG C 300-OT progenitor stars 
dPrieto et alj 120081 iPrietol 120081) . This is simply a coinci- 
dence driven by the present luminosity and dust radius, be- 
cause in steady state a surviving star could not support the 
necessary optical depths given the velocities needed to have 
new material at such large d istan ces. The populati on stud- 
ies bv IThompson et all d2009l) and lKhan et all d2010l) imply a 
lifetime in the obscured phase of order t B ~ 10 4 years. The 
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while the star is in its obscured phase. No star could sus- 
tain this for v w ~ 10 3 km/s, so the current state cannot rep- 
resent a long lived period and must be a transient phase. 
The progenitors , unde r the super- AGB star hypothesis of 
Thompson et alj (120091) . can have the necessary lifetimes be- 
cause dust driven wind velociti es of v w ~ 10-20 km/s are so 
low (e.g. Ilvezic & ElitzurJl2010h . However, such a slow wind 
commencing after the transient would not even have started to 
form dust at this point in time. 

The second possibility is the traditional view of the super- 
nova impostors as stars that briefly enter a high mass loss 
state during the transient to produce a shell of ejected materia l 
that then forms dust. As discussed in lKochanek et alJ (12012b . 
SN 2002bu probably started forming dust too soon after the 
transient peak to have formed it in ejected material. The more 
serious problem is that between 2004 and 2012 an expanding, 
dense shell should have a radius growing as Rj„ oc t and an op- 
tical depth dropping as Ty oc ?~ 2 . If we take the elapsed time 
from discovery to the IRAC observations in 2004 (736 days) 
and 2012 (3566 days), the radius should have increased by a 
factor of 4.9 and the optical depth should have dropped by a 
factor of 23.5. Fig. [6] illustrates this by comparing the mean 
expansion velocities and optical depths for the two epochs 
and both grain types. If we scale the values in 2004 to their 
expected values in 2012 simply using the expected temporal 



5 




Wavelength [/u.m] 



FIG. 4. — The top left panel shows the SN 2002bu SED roughly 2 years after the explosion, while the other 4 panels show the current SED for sources A, B, 
and C (see Fig. [2), in the case of source A also treating HST near-IR measurements as upper limits (top right panel). Magnitudes from Table[T]are converted to 
fluxes, and then to luminosities as L = 4ttD 2 uF 1 , where D = 5.8 Mpc. The SED models are black bodies with Ay = 0.012 mag of total (Galactic) extinction. The 
thick black curves show the probability averaged SEDs for graphitic (solid) and silicate (dashed) dust, while the th in curves show the spr ead around the mean 
values. In the case of the 2004 SED, solid and dashed lines are the best fit graphitic and silicate models taken from Kochanek et al. 12012) and correspond to a 
T* = 20000 K, L* ~ 10 5 ' 9 Lq star surrounded by an optically thick (tv ~ 30) shell of dust expanding at ~900 km/s. 

to a dense pre-existing wind. In this scenario, the progen- 
itors were obscured by a dense, dusty wind where the dust 
is destroyed by the shock breakout luminosity from the ex- 
plosive transient so as to leave the transient with little obscu- 
ration at peak. The densities are so high, however, that the 
dust can then reform and re-obscure the system, consistent 
with the early dust formation suggested by the spectra and 
photometry. In the SN 2008S scenario, the material already 
exists at these distances and so dust can begin (re-)forming 
very rapidly. If the luminosity is then powered by a combina- 
tion of any surviving star and an expanding shock wave, the 
optical depth to the shock front is only dropping as ry oc l/t . 
The resulting drop in the optical depth by only a factor of 
4.9 is far more compatible with the observations, as shown 
in Fig. [6] Arguably, the shock should also be slowing as it 
expands through the dense wind. Self-similar solutions (e.g. 
I Chevalier! (11983b ') give expansion rates of R oc f and v s oc 
where s = (n - 3)/(n - 2) for ejecta with an effective density 
profile p oc R~" expanding into a p oc R~ 2 wind. The typical ap- 
proximations for the structure of the ejecta are n = 7 (s = 0.8) 
and « = 12 (s = 0. 1), leading to a modest slowing with time. In 
Fig. [6] we show the scaling of the conditions in 2004 to 2012 
for this shock scenario (v s oc f~ 0,2 , R oc t °' 8 , r oc t ~ ' 8 ). 

If a significant part of the present day luminosity is driven 
by X-rays from an expanding shock, inferences about the 
temperature of the illuminating source and the exact value 
of the optical depth become somewhat problematic. In the 
shock scenario, the X-rays are absorbed by the dense gas and 
remitted as a complex, non-thermal (mainly) emission line 
spectrum, which is then absorbed and reradiated by the d ust. 
While the shock emission models of lAllen et aD (120081) do 
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FIG. 5.— The mid-IR CMD of NGC 2403 (Khan et al. 2010). The color 
evolution of SN 2002bu is marked with black squares connected by a dotted 
line, where the top point corresponds to the 2004 epoch, while the next to 
2011 and 2012 epochs (see Table ff}. For comparison, we show the color 
evolution of SN 2008S from 2008 to 2012 (blue circ les) and NGC 300 -OT 
from 2008 to 2012 (red triangles) as summarized in Kochanek (2 0111) and 
ISzczygiel et ST] <2TjT2l) . 



scalings, we see that the observed properties in 2012 have a 
significantly higher optical depth and may require some de- 
celeration of the expansion rate. 

This leaves the scenario developed by Kochanek] d201 ll) for 
SN 2008S and the NGC 300-OT, where the absorption is due 
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FIG. 6. — The average expansion velocity and optical depth of the DUSTY models for SN 2002bu for graphitic (left) and silicate (right) dusts. The black (red) 
contours encompass the la region for two parameters for the 2004 (2012) epoch. The dashed green contours (labeled shock) show the expected parameters in 
2012 given the parameters in 2004 assuming the radius is expanding as R oc t" with x = 0.8 (left) or 1.0 (right) and the optical depth is dropping as ry oc \/R. 
The dotted blue (labeled shell) shows the expected parameters in 2012 if the radius is expanding as R oc t and the optical depth is dropping as Ty oc 1/r 2 . Shock 
models to the right of the heavy solid line are excluded (Icr) by the upper limits on the X-ray flux in the Swift observations. If the HST sources are treated as 
upper limits, the 2012 contours expand upwards to higher optical depths. The optical depth is related to the wind density by Eqn.[T]where we used m = 1.7 (0.84) 
for graphitic (silicate) dusts. 



not extend to the regimes considered here, they generally have 
relatively low near-IR emission compared to optical/UV emis- 
sion and so would resemble a hot star model. The estimates of 
the optical depth, which are driven by the near-IR detections, 
should also be viewed as being only logarithmically correct 
in this scenario. The overall luminosity, dust temperature and 
dust radius, which are largely determined by the mid-IR emis- 
sions from the dust, should still be accurate. 

If we fit the expansion of the dust radius as a power law, 
R oc t x , we find x = 0.67 (0.26 < x < 1 .05, la) for the graphitic 
models and x = 0.54 (0.29 < x < 0.80, la) for the silicate 
models. The silicate models are consistent with linear expan- 
sion only at 2a. Similarly, if we fit the evolution of the optical 
depth as a power law, Ty oc t y , we find y = -1.27 (-1.69 < 
y < -0.90, la) and y = -0.79 (-1.29 < y < -0.47, lcr) for 
the graphitic and silicate models, respectively. The Ty oc 1 jt 2 
scaling of expanding shells is ruled out at roughly 2a. We 
can also fit the evolution of the optical depth as a power law 
in radius, Ty oc R z , finding z = -1.75 (-3.90 < z < -1.01) and 
z = -1.52 (-3.46 < z < -0.69) for the two models. While 
this appears to be consistent with the Ty oc R~ 2 scaling for an 
expanding shell, there is a strong covariance between this ex- 
ponent and the expansion rate of the shell - a solution with 
Ty oc R~ 2 must also roughly have R oc t l l 2 . Essentially, for the 
same change in optical depth there is considerable uncertainty 
in the radius because we lack a complete dust SED in 2012 
for constraining the dust temperature. In the two-dimensional 
probability distribution of the x and z exponents, the expand- 
ing shell solution with x = 1 and z = -2 is still ruled out at 
roughly 2cr. The evolution is consistent with the Ty oc l/R 
and R oc f° 8-1 evolution of the expanding shock model. 

A shock moving through a wind at velocity v s has a charac- 



teristic X-ray energy of 
16 



E s = ^nipV 2 = 1.2 



1000 km/s 



keV 



(3) 



where the mean molecular weight is /i = 0.6, and produces 
luminosity 

L S = \U (4) 

2 v w 

where e < 1 is the efficienc y with wh ich the shock energy 
is radiated as X-rays (e.g. IChevalien 11981 [Chugail 1 1 9921 
Chuga i & Danzigerl Tl994, Ch evalier & Franssonlll994l) . The 
radius of the shock R = v s t = 3.1 x lO'^v./lOOO km/s) cm 
is simply related to its velocity, where t = 3566 days is the 
elapsed time from the transient to the X-ray observation. This 
means that the hydrogen column density outside the shock ra- 
dius is 



N H 



M 



4TTV w R(lAm p ) 



= 7.5 x 10 



21 



M 



10- 4 M Q /year 



10 km/s \ / 1000 km/s \ 
v s J 



(5) 



cm 



Combining these three equations with the PIMMS models for 
absorbed thermal bremsstrahlung, we can determine the range 
of shock velocities and wind density parameters (M/v w ) that 
would violate the Swift limits, as shown in Fig. [6] If the shock 
velocity is as high as the v s = 893 km/s FWHM reported by 
iSmith et ail ( 1201 11) . then there are allowed solutions for low 
wind densities (where the X-ray luminosity is low but little 
absorbed), and high wind densities (where the X-ray lumi- 
nosity is high but heavily absorbed). If the shock velocity 
v s < 600 km/s, then the X-ray emission limits are satisfied for 
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any wind density. The parameters derived from the photo- 
metric fits broadly satisfy these limits on the X-ray emission. 
Half the graphitic solutions violate the limit in Fig. [6] but the 
limit as drawn is for 100% conversion of the shock energy into 
soft X-rays and a la detection threshold, which is a rather op- 
timistic representation of the detection threshold. The basic 
picture of Fig.[6]is little changed if we use the larger distances 
to NGC 4242, with all solutions and limits simply shifting to 
higher velocities. While this increases the shock luminosity, 
the net effect is limited because the X-ray luminosity limits 
from the observations become correspondingly weaker. 

4. SUMMARY 

Combining Spitzer, Hubble and Swift observations of 
SN 20 02b u, we confirm th e arguments in Thompson et al. 
(2009) and Koc hanek et al.l d2012l) that this source is a mem- 
ber of the SN 2008S class of transients. While we have no 
direct evidence that the progenitor was self-obscured, dust ap- 
pears to have begun forming before any ejecta would be cool 
enough to do so and the optical depth of the dust is dropping 
significantly more slowly than the -ry oc 1 jt 2 scaling of an ex- 
panding shell of material. The evoluti on is broadly con sis- 
tent with the expanding shock model of Kochanek] d201 lb for 
SN 2008S and the NGC 300-OT. It is not consistent with an 
expanding shell of ejected material or a reformed, slow wind 
around a surviving dusty star. Unfortunately, our Swift X-ray 
observations only constrain the allowed parameter space since 
w e failed to detec t X-ray emission from the shock. 

iGogarten et al.l (120091) argued from an analysis of the stellar 
populations near the NGC 300-O T that the progenito r could 
be as massive as ~ 2OM , while Prieto et al. (2008) argued 
for ~ 10M© because the progenitor appeared to be an extreme 
AGB star based on its location in mid-IR CMDs. The ambi- 
guity is possible because the most luminous AGB stars have 
luminosities comparable to evolved ~ 20M Q stars. Formally, 
however, these population analyses only provide upper mass 
limits, so it is entirely possible to have a transient (even a su- 
pernova) from a ~ 1OM star in a region containing ~ 20M Q 
stars. However, the lo cal stellar environmen t of SN 2002bu, 
like that of SN 2008S dSzczygiel et al.ll2012l) . appears to con- 
tain no stars more massive than ~ 10M Q . In fact, the envi- 
ronment of SN 2002bu appears to favor still lower masses, 
closer to 5M Q than 10M Q . This strongly favors the inter- 



pretation of lPrietoetail d2008l) and iThompson etaTI (120091) 
that the SN 2008S class of transients is associated with AGB 
stars. The very low mass scale implied by the environment of 
SN 2002bu does not, however, favor the po ssibility that the 
transients are electron capture supernovae (Thompso n et al.l 
120091 iBotticell a et alj|2009|). as these are only expected for 
masses of 8-10M Q (e.g. iPoelarends et al.ll2008l) . Obtaining 
deep optical images of the SN 2008S and 2002bu fields to 
better constrain the extinction and temperatures of the nearby 
stars would greatly improve the characterization of the stellar 
populations over our present near-IR results. 

As with all the "supernova impostors", the only way to de- 
termine their natures is to continue to track the evolution of 
th eir spectral energy dist ributions. It is clear from our survey 
in lKoch anek et all (l2012l) that there are at least two classes of 
sources and that no sources with adequate data were consis- 
tent with a brief period of high mass loss during the transient 
leading to the formation of a dusty expanding shell that ob- 
scured the source. The strength of this conclusion is strongly 
limited by the fragmentary nature of the data in both wave- 
length and time. Here we see that the SN 2008S transient 
class can remain dust-obscured for long periods of time, pri- 
marily because the optical depth is not evolving like that of an 
expanding shell. However, if object A is the near-IR counter- 
part of SN 2002bu, the veil is beginning to thin. 
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TABLE 1 
Photometry of SN 2002bu 



Date (UT) 


MJD 


Pi/Program 


F435W 


F555W 


F625W 


F814W 


FUOW 


F16QW 


[3.6] 


[4.5] 


[5.8] 


[8.0] 


[24.0] 


[70.0] 


Comments 








[mag] 


[mag] 


[mag] 


[mag] 


[mag] 


[mag] 


[jUJy] 


[/Jy] 


[MJy] 




[/«Jy] 


[/Jy] 




2004-04-25 


53120 


Fazio/69 






















750 ±26 


<4200 


1 


2004-05-02 


53128 


Fazio/69 














544 ±7 


806 ±9 


1098 ±12 


1204 ±20 






1 


2005-03-20 
2011-06-23 


53456 
55735 


Filippenko/10272 
Kochanek/80015 


< 25.67 


< 25.25 


< 24.83 


< 24.40 






8±3 


11±1 










1 


2012-01-31 


55958 


Kochanek/80015 














5±1 


13± 1 












2012-02-15 


55972 


Kochanek/12450 










25.64 ±0.10 


24.85 ±0.13 














A 



27.21 ±0.39 26.72±0.71 B 
25.78±0.12 25.53±0.26 C 

NOTE. — All the magnitude upper limits are 3o\ The date of discovery of the transient is MJD 52392.3. (1) From Kochanek et al. 1 2012); (A) measurement for Source A; (B) measurement for Source B; (C) 
measurement for Source C. 
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TABLE 2 

Spectral Energy Distribution Models 



MJD 


iog 10 r» 

(K) 


log 10 L« 

(£©) 


(K) 


log TV 


log 10 R,„ 
(cm) 


log io v > 
(km/s) 


Comment 


53128 
53128 
55972 
55972 


3.99 ±0.31 
3.88 ±0.32 
4.32 ±0.09 
4.43 ±0.04 


4.92 ±0.01 
4.95 ±0.02 
4.52 ±0.22 
4.95±0.18 


3.04 ±0.07 
3.09 ±0.10 
2.65 ±0.05 
2.77 ±0.04 


1.58 ±0.21 
1.77 ±0.27 
0.70 ±0.25 
1.21 ±0.11 


15.77±0.13 
15.62 ±0.07 
16.18 ±0.25 
15.98±0.17 


2.94 ±0.1 3 
2.80 ±0.07 
2.72 ±0.25 
2.52 ±0.17 


Graphitic 
Silicate 
Graphitic 
Silicate 



NOTE. — These include a weak prior on the stellar temperature, log 10 T* = 4.0 ± 0.3 (3000 < T* < 30000), and the velocity 
log 10 (v s /km/s) = log 10 893 ± 0.3 at which R m expands. These models all assume that source A is the near-IR counterpart of 
SN 2002bu. If it is a chance coincidence, then the optical depth scale can be significantly higher. 



